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ABSTRACT 

Stimulated  by  recent  experimental  observations  of  room  temperature  ferromagnetism  of 
Mn^-Cdi-^GePa,  we  investigate  the  structural,  electronic  and  magnetic  properties  of  these 
systems  as  a function  of  Mn  concentration  by  means  of  first-principles  donsity-functional- 
thcory-based  codes.  Moreover,  we  investigate  the  effect  of  the  anion  substitution  (P  vs  As) 
in  Mn-rich  chalcopyirites.  Our  calculations  indicate  that  the  antiferromagnetic  alignment  is 
the  most  stable  ordering  for  all  the  systems  studied,  at  variance  with  that  experimentally 
reported.  Moreover,  we  find  that  there  is  a slight  reduction  of  the  total  magnetic  moment, 
per  Mn  atom  from  ~5  pn  in  all  the  Cd-rich  P-based  chalcopyrites  to  ~4  fin  in  the  Mil  rich 
MnGeP2  and  MnGeAsa  systems. 

INTRODUCTION  AND  COMPUTATIONAL  DETAILS 

A great  effort  is  being  devoted  to  explore  new,  advanced  magnetic  semiconductors  for 
spintronics  applicat,ions[l],  based  on  crystal  structures  that  transcend  the  limitations  (defect 
formation  and  too  low  operating  temperatures)  of  the  magnetic  zinc-blende  systems  explored 
so  far  ( e.g . MiijGai-^As).  Within  this  framework,  the  chalcopyrites  seem  to  be  good 
candidates;  the  expected  advantage  of  these  systems  is  that  Mn  can  readily  substitute  for 
the  typo-II  cations,  as  has  been  demonstrated  for  IR^Mn^VI  alloys  with  x up  to  1.0  without, 
defect  formation,  due  to  the  natural  tendency  of  Mn  to  adopt  a +2  valence  state. 

Recently,  Medvedkin  et  al.  [2]  deposited  a.  Mn  overlayer  in  vacuum  on  the  surface  region 
of  a CdGeP2  crystal,  followed  by  a solid  phase  reaction  at  elevated  temperatures.  A rele- 
vant finding  of  Medvedkin  et  al.  is  the  room  temperature  ferromagnetism  in  highly  doped 
Mn^Cdi _IGeP2  - which  constitutes  a tremendous  improvement  from  the  Tc;  of  110  K found 
in  Mn^Gai-jAs  at  x = 5.3  %.  Moreover,  their  photoluminescence  spectrum  at  20  K shows 
a peak  around  3.2  eV,  from  which  they  proposed  a Mn  induced  enhancement  of  the  energy 
gap  of  CdGeP2  (Eg  = 1.83  eV).  This  is  quite  different,  from  the  III-V  alloys,  in  which  Mn  is 
known  to  reduce  the  band-gap  [7]. 

In  order  to  understand  the  magnetic  interaction  in  chalcopyritc  semiconductors,  we  inves- 
tigated MnICdi_«GeP2  for  various  x and  MnGo(V)2  whore  V = As,  P with  the  full-potential 
linearized  augmented  plane  wave  (FLAPW)[3]  method.  This  is  one  of  the  most  accurate  ab- 
initio  methods,  in  which  there  is  no  artificial  shape  approximation  for  the  wave  functions, 
charge  density  and  potential.  For  all  atoms,  the  core  and  valence  states  are  treated  fully— 
and  seini  relativistically  ( i.e . without  spin  orbit  coupling),  respectively.  The  muffin  tin 
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(MT)  radii  for  Cd,  Ge,  Mn  are  set,  to  2.3  a.u.;  2.0  and  1.8  a.u.  arc  used  for  As  and  P, 
respectively.  An  energy  cut  off  of  9.0  Ry  was  employed  for  the  wave  function  expansion  in 
the  interstitial  region,  whereas  a 49  Ry  cut-off  was  used  for  the  charge  density  and  potential. 

For  all  the  systems  studied,  we  used  the  generalized  gradient  approximation  (GGA) 
functional  of  Perdew-Burke-Ernzerhof  [5]  rather  than  the  bare  local  density  approximation 
(LIDA)  exchange-correlation  functional  [G],  within  density  functional  theory;  the  GGA  was 
proved  to  give  better  agreement  with  available  experimental  structural  data  for  Mn  pnictides 
[8]  and  is  therefore  likely  to  describe  more  accurately  the  Mn  rich  chalcopyrites. 

STRUCTURAL  DETAILS 

Let  us  first  focus  on  the  structural  details.  Recall  that  the  chalcopyrite  differs  from  the 
zincblende  crystal  structure  by  a doubling  of  the  unit  cell  along  a four  fold  axis,  rendering 
the  system  body  centered  tetragonal.  Through  total  energy  calculations,  we  fully  optimized 
the  MnGeAs2  system,  obtaining  u = 5.83  A,  c/a  = 1.95±0.05  and  u ~ 0.25.  The  estimated 
numerical  uncertainty  on  the  ?/  value  can  be  ascribed  to  the  very  small  effect  on  the  total 
energy  due  to  the  tetragonal  deformation  (i.e.  7/  = c/2a  ^ 1)  around  the  ;/  equilibrium 
value.  Our  optimized  structural  values  are  in  good  agreement  (within  0.5  % on  the  a lattice 
constant.)  with  tin;  values  predicted  by  the  “conservation  of  tetrahedral  bonds”  (CTB)  plus 
?/  = »/ irf”  theory  discussed  in  Refs.  [9].  In  our  work,  Pauling’s  tetrahedral  radii  were  chosen 
as  1.48,  1.22,  1.22  and  1.10  A for  Cd,  Ge,  As  and  P respectively;  the  Mn  tetrahedral  radius, 
derived  from  the  Mn  chalcogenides  is  1.41  A.  With  these  values,  we  obtained  a = 5.86  A 
and  c./a  ~ 1.9.  The  good  agreement  for  the  As-based  chalcopyrite  suggests  that  the  “CTB 
plus  7/  = 7/(rt”  rule  may  even  hold  for  the  similar  P-based  MnGeP2  system;  we  therefore 
used  a = 5.67  A,  c/a  = 1.9.  Regarding  the  MnCdGeP2  structures[10],  we  interpolated  the 
lattice  constants  (a  and  if)  from  the  experimental  values  for  CdGeP2  and  the  ”CTB  plus 
V = 7/irt”  value  for  MnGeP2.  This  approximation  is  equivalent  to  Vegard’s  law  and  works 
well  in  this  system.  For  instance,  it  is  reported  by  Medvedkin  et  al.  [2]  that  the  lattice 
constant  decreases  to  5.7  A when  the  Mn  concentration  becomes  greater  than  Cd,  which 
means  x in  the  range  of  50  % to  53  % since  the  maximum  x is  53  % in  the  experiment.  Our 
interpolated  latt  ice  constant  for  a:=0.5  is  5.71,  which  is  different  from  experiment  by  only  0.4 
%.  Moreover,  in  the  MnCdGeP2  systems,  we  optimized  the  internal  degrees  of  freedom  using 
both  the  I)mol3[4]  ( i.c . density  functional  theory  for  molecules  and  3-dimensional  periodic 
solids)  and  the  FLAPW  codes.  More  accurate  details  about  the  structural  properties  of 
Mn:CdGeP2  can  be  found  in  Ref.  [10]. 

We  focused  on  the  dependence  of  the  electronic  and  magnetic  properties  in  semiconduct- 
ing chalcopyrites  as  a function  of  i)  the  type  V anion  (i.e.  P vs  As)  in  Mil  rich  systems  and 
ii)  the  Mn  concentration  (x  = 0.25,  0.5  and  1). 

THE  MnGcV2  (V  = As,  P)  SYSTEMS 

Let  us  first,  discuss  the  effects  due  to  the  anion  species  by  comparing  some  relevant 
properties  in  MnGeP2  and  MnGeAs2  shown  in  Table  1.  As  for  the  antiferromagnetic  (AFM) 
configurations,  we  considered  a 1x1  chalcopyrite  “spin  superlatt.ie.e”  along  the  [001]  direction. 
We  find  that  both  the  P-  and  As-based  systems  largely  prefer  the  AFM  alignment,  since  the 
AFM  is  mill'll  lower  in  energy  than  the  ferromagnetic  (FM)  structures.  Let  us  now  consider 
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Table  1:  Comparison  of  relevant  structural  and  magnetic  properties  of  MnGeP2  and 
MnGeAs2. 


MnGeP2 

MnGeA.s2 

FM 

AFM 

FM 

AFM 

Total  energy  (meV/Mn) 

0 

-160 

0 

-60 

/i Mr,  (Bohr  magn.) 

3.55 

3.55 

3.74 

3.68 

Total  spin  moment/Mn  (Bohr  magn.) 

4.0 

- 

4.3 

- 

Mn-anion  bond  length  (A) 

2.44 

2.43 

2.52 

2.51 

Figure  1:  Spin  density  of  ferromagnetic  MnGcAs2,  plotted  in  the  [110]  plane.  Positive 
(negative)  spin-density  is  represented  by  solid  (dashed)  lines. 


Energy  (eV) 

Figure  2:  Total  density  of  states  of  MnGeP2  (dashed  line)  and  MnGeAs2  (solid  line)  for  FM 
(upper  panel)  and  AFM  (lower  panel)  ordering.  The  zero  of  the  energy  scale  denotes  the 
Fermi  level. 

the  magnetic  moments.  In  this  case,  Mil  substitutes  a type-II  cation;  we  would  therefore 
expect  a magnetic  moment  close  to  5.0  pD.  However,  note  that  the  local  moment  relative  to 
the  free?  Mil  atom  (/«o  = 5 /iB  in  the  rf5  configuration)  is  reduced  due  to  the  mixing  of  the 
unoccupied  Mn  rf*  bands  with  the  occupied  p bands  by  the  p — d hybridization.  Note  that  in 
both  As-  and  P-based  cJialcopyrit.cs,  the  Mn  induced  FM  moment  at  Ge  sites  gives  positive 
contributions,  while  both  P and  As  have  a very  small  induced  magnetic  moment  within  the 
muffin  tin  sphere.  As  an  example,  we  plot  in  Fig.  1 the  spin-density  in  MnGeAs2  in  the 
[110]  plane;  it  is  evident  that  a negative  spin  density  is  present  around  the  As  atom,  even 
though  there  is  a positive  spin  density  along  the  bond. 

Let  us  now  discuss  the  electronic  properties  of  the  MnGeAs2  and  MnGeP2  systems.  In 
Fig.  2 the  spin-resolved  total  density  of  states  for  the  FM  and  AFM  alignments  is  reported. 
Let  us  first  focus  on  the  FM  alignment,  whose  main  features  arc:  i)  around  -12  eV  we  find 
the  hybridized  Ge  s and,  mostly,  anion  s bands;  ii)  the  structure  around  -7  cV  has  a similar 


— - MnGeP2 
MnGeAs2 
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Table  2:  Relevant  structural,  electronic  and  magnetic  properties  in  Mnx.Cdi_xGeP2  for  var- 
ious Mn  concentrations.  The  x = 0.5  and  1.0  are  lxl  ” spin-superlattices”  ordered  along  the 
[001]  direction,  whereas  the  x = 0.25  is  ordered  along  the  [100]  direction. 


X 

[ 0.25 

0.5 

1.0 

FM 

AFM 

FM 

AFM 

FM 

AFM 

MnP  bond  length  (A) 

2.43 

2.43 

2.44 

2.43 

2.44 

2.43 

Total  energy  (meV /Mn) 

0 

-30 

0 

-35 

0 

-160 

Total  spin  moment/Mn  (Bohr  magn.) 

- 

~5 

- 

- 

Energy  gap  (eV) 

0.1C 

0.58 

- 

0.49 

- 

0.71 

origin,  with  a more  marked  Go  s character;  Hi)  the  states  at  higher  energy  (i.e  from  -5 
eV  to  +5  oV)  have  a mixed  character.  In  particular,  it  is  evident  that  the  Mn  $ states 
are  almost  totally  occupied  (c.f.  feature  peaked  at  around  -3.4  eV  in  the  spin  up  DOS), 
whereas  the  d1,  states  are  mainly  unoccupied  and  hybridized  with  the  p anion  states,  (c.f. 
feature  going  from  -0.8  eV,  peaked  at  around  0.5  eV  and  then  up  to  5 eV  in  the  spin  down 
DOS).  The  pronounced  feature  in  the  spin-  down  DOS  from  -5  eV  to  -0.8  eV  is  mainly  due 
to  hybridization  between  Ge  p and  anion  p states.  We  can  estimate  a d exchange  splitting 
of  about  4 cV.  Neither  the  spin-up  nor  the  spin-down  density  of  states  show  a gap,  therefore 
resulting  in  a metallic  system;  on  the  other  hand,  there  is  a small  gap  in  the  spin-down 
DOS  at  around  -0.8  eV  (between  the  Gc-anion  p bands  and  the  Mn  d states)  which  could 
lead  to  half-met.allicity,  possibly  by  proper  modifications  in  the  lattice  constants.  The  main 
difference  between  As  and  P anions  stems  from  the  higher  hybridization  [8]  with  the  Mn  d 
bands,  shown  by  an  almost  rigid  shift  of  the  P (Mn)  features  towards  lower  (higher)  binding 
energies  with  respect  to  the  As  case.  As  for  the  AFM  DOS,  we  can  outline  features  similar  to 
the  FM  case.  We  still  find  a metallic  character,  even  though  a depletion  of  states  is  present 
just  around  the  Fermi  level,  suggesting  a starting  gap  opening  effect.  This  difference  with 
respect  to  the  FM  case  is  related  to  the  i)  less  evident  hybridization  between  Mn  eft  with 
anion  p'1  states  and  ii)  totally  unoccupied  Mn  # antibonding  states. 

THE  Mn3,Cd1_xGeP2  STRUCTURES 

Let  us  now  focus  on  the  effect  of  Mn  vs  Cd  substitution  for  Cd  in  the  P based  chalcopy- 
rites.  We  report  in  Table  2 the  relevant  structural,  electronic  and  magnetic  properties  of 
these  systems.  First  of  all,  note  that  the.  AFM  alignment  is  the  most  stable  for  all  Mn  con- 
centrations, although  the  energy  difference,  between  FM  and  AFM  configurations  is  strongly 
reduced  with  respect  to  the  Mn-rich  chalcopyrite.  These  results  are  not  in  agreement  with 
the  experimental  observations[10],  that  found  MnxCd|_xGeP2  to  be  a room-temperature  fer- 
romagnet.  The  reason  for  this  disagreement  is  not  clear.  In  our  calculations,  Mn.,.Cd]_xGeP2 
is  an  ideal  periodic  bulk  structure,  while  in  experiment  it  is  a surface  doped  system  with 
Mn  concentration  that  decays  from  the  surface  to  the  bulk.  To  resolve  this  discrepancy, 
a more  thorough  and  careful  experimental  characterization  is  called  for.  As  for  the  total 
spin  moment,  we  note  that  the  Mn  induced  moment  on  the  Cd  atom  is  positive  and  this  is 
likely  to  justify  the  higher  total  moment  for  x = 0.25,  0.5  than  for  x = 1.0  systems. 


T6.9.5 


Finally,  we  focus  on  the  energy  gap  in  Mii:rCdi_ICGeP2.  Recall  that  within  GGA  we  obtain 
strongly  underestimated  band  gaps,  with  respect  to  experimental  values;  for  example,  we 
obtain  0.93  eV  for  the  CdGeP2  crystal,  to  be  compared  with  the  experimental  value  of  1.83 
eV.  In  the  AFM  state,  the  energy  gap  decreases  to  0.49  eV  for  x — 0.5,  although  it  goes  up 
to  0.71  eV  for  x — 1.0.  In  all  cases,  the  energy  gap  of  Mn:CdGeP2  does  not  exceed  that  of 
pure  CdGeP2.  In  the  FM  state,  there  is  no  energy  gap  except  for  x = 0.25.  Therefore,  the 
effect  of  Mil  doping  on  the  energy  gap  in  CdGeP2  is  similar  to  that  in  Mn^-Ga^  As,  where 
a reduction  of  the  energy  gap  of  the  host  semiconductor  is  found.  This  disagrees  with  the 
experimental  observation  of  Medvedkin  et  «/.[2],  showing  an  increase  of  the  energy  gap  upon 
incorporation  in  CdGeP2. 

CONCLUSIONS 

The  recently  reported  room  temperature  ferromagnetism  in  Cdi_;,,Mn,.GeP2  was  investi- 
gated for  x — 1.0,  0.5  and  0.25  using  the  FLAPW  and  DMol  codes,  within  the  generalized 
gradient  approximation  to  the  density  functional  theory.  Moreover,  we  have  compared  the 
MnGeA&j  and  MnGeP2  systems,  in  terms  of  structural,  electronic  and  magnetic  properties. 
At  variance  with  recent,  experimental  observations,  we  find  that  the  total  energy  of  the  AFM 
state  is  always  lower  than  the  corresponding  FM  state  for  all  the  systems  studied.  The  total 
magnetic  moment  per  Mn  atom  is  close  to  5 //y  for  the  Cd-rich  P-based  systems,  whereas 
it  is  reduced  in  the  Mn-rich  P-  and  As-based  chalcojryrit.es. 
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